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Surface Spectral Function of Momentum-dependent Pairing Potentials in a 
Topological Insulator: Application to Cua.Bi2Se3 

Liang Chen and Shaolong Wan* 

Institute for Theoretical Physics and Department of Modern Physics 
University of Science and Technology of China, Hefei, 230026, People's Republic of China 

(Dated: February 26, 2013) 

We propose three possible momentum-dependent pairing potentials for candidate of topological 
superconductor (for example Cua;Bi2Se3), and calculate the surface spectral function and surface 
density of state with these pairing potentials. We find that the first two can give the same spectral 
functions as the fully-gapped and node-contacted pairing potentials given in [Phys. Rev. Lett. 105, 
097001], and that the third one can obtain topological non-trivial case which exists fiat Andreev 
bound state and preserves the C3 rotation symmetry. We hope our proposals and results be judged 
by future experiment. 
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Recently, topological insulators [1-3] have attracted 
great attention in condensed matter physics for their 
physical properties. The experimental [4-!)] and theo- 
retical [10-14] investigations show that topological insu- 
lators (TI) are fully gapped in the bulk but gapless on 
the surface, which is protected by time-reversal symme- 
try. And these surface states indicate the massless Dirac 
fermions. Then, the researches on topological insulators 
are generalized to on topological superconductors (TSC) 
[15-17] soon. Similarly, researches show that topological 
superconductors also have gapless surface states, which 
indicate massless Majorana fermions [18], and that the 
property on TSC is protected by topological bulk proper- 
ties and characterized by a topological invariant [L5, 19]. 
As we know, the Majorana fermions are of great interest 
in fundamental physics and have potential applications 
in quantum computation [20, 2J]. 

The experiment [22] finds that a superconductive phase 
is induced at transition temperature Tc ~ 3.8A' when 
copper atoms are doped into topological insulator Bi2Se3 
with the concentration of Cu in range 0.12~0.15. And 
the work [23] shows that the surface state of CuxBi2Se3 
is topological non-trivial. A recent experiment [24] fur- 
ther confirms the existence of topological surface state 
by measuring the surface density of states (SDOS). How- 
ever, a more detailed analysis show that a gapless and 
topological non-trivial bulk band structure may be pre- 
ferred. 

Up to now, the exact pairing mechanism still be un- 
clear, but there are some theoretical proposals on the 
pairing symmetry of Cu2;Bi2Se3, including odd-parity 
pairing potential proposed by Fu and Berg [25] and Sato 
[26, 27] and the explanation of Sasaki et.al.\lA\ accord- 
ing to experimental results. However, in all theoreti- 



cal proposals, they assume that the pairing potentials 
are momentum-independent. For this reason, we ask 
whether we can search for some pairing symmetries which 
induce new topological surface state and act as candi- 
dates of pairing symmetry of Cu2;Bi2Se3, if the pairing 
potentials are momentum-dependent. Further more, we 
know that the A4 suggested by Ref. [24] breaks C3 ro- 
tation symmetry of the rhombohedral lattice. So we also 
want to know whether there is possible pairing poten- 
tial which is topological non-trivial, node-contacted and 
preserves the C3 rotation symmetry. In order to search 
for answers of these questions, in this paper, we pro- 
pose three momentum-dependent pairing potentials for 
Cuj;Bi2Se3, and calculate the surface spectral functions 
with these pairing potentials. And we find that the first 
two sorts of our pairing potentials are similar to A2 and 
A4 proposed in Ref. [25] and can obtain what A2 and 
A4 give, and that the third one can get topological non- 
trivial case which exists fiat Andreev bound state (ABS) 
and preserve the C3 rotation symmetry. 

As reported [23], near the F-point, the band dispersion 
of normal state of Cu2,Bi2Se3 can be described by the 
low energy effective Hamiltonian [14] for Bi2Se3, with a 
finite chemical potential in the conduction band induced 
by copper doping. The Hamiltonian is 

h{k) = {Mt, -h)+t^ (A(fc^CT^ -I- kyay) + Bk,a,) , (1) 

where M is the rest mass, n is the chemical potential, 
A and B are Fermi velocity along different directions, 
T2 = ±1 denotes the two orbits, and crx,y,z are spin Pauli 
matrices. For the superconductivity phase, the Hamilto- 
nian can be written in the Bogoliubov-de Gennes (BdG) 
formalism: 



n 



h{k) 

At(fe) 



A(fe) 
-K^{-k) 



(2) 
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where A(fe) is the pairing potential (4x4 matrix). For the 
time reversal invariant case, the pairing potential can be 



divided into two parts according to inversion symmetry: 



A(fe) = Ai(fc) + A2(fe), 



Ai(fc) = -ifcjCTy (cTaAi-'g + r^(TQAi_-'^ 



A2(fe) 



k. 



[Tx(Ja 



^2 



ir^A^^ 



(3) 
(4) 

(5) 



where A"-'^, A"-'^^, Aj"* and A'^'' are real functions 
of momentum and inversion symmetric, the summa- 
tion convention is used in this paper. One can find 
that A2(fc) = V~^ A2{k)V is inversion symmetric and 
Ai(fc) = —V^^Ai{k)V is inversion anti-symmetric, 
where the inversion operator is 7-' = Tz in coordinate 
space. As we know, if the pairing potential is dominated 
by Ai(/c) and fully gapped in the Brillouin zone, the cri- 
terion for topological odd-parity superconductor [25-27] 
claims that the system is topological nontrivial. 

In the following, we consider three typical cases of 
Ai(fc) and calculate the surface spectral function with 
them. In order to calculate the spectral function numer- 
ically, we consider a lattice model which the low energy 
effective Hamiltonian is Eq.(l). For the normal state of 
the Hamiltonian, we use the model and parameters given 
in supplemental material of Ref. [24] . 

In the first case, for superconductive potential, Eqs.(3)- 
(5), we consider A^*^^^^ = A^^ = A^^ = and Ai,, = 
Adiag(yl, ~A, B), as an example, and other cases are sim- 
ilar. In this case, the pairing potential takes 



A(fc) = A 






Bk^ 



®ro, (6) 



where A is a dimensionless parameter determined by en- 
ergy gap of superconductivity, tq is a 2 x 2 identity matrix 
in orbital space. One of compactifications of Eq.(6) can 
be given as 



A(A;) = A 



-A, Ai 



Ai 



At 



®To, 



(7) 



where we refer to the definition of ^i, A^ in the supple- 
mental Material of Ref. [24] , and take the same param- 
eter values of Ref. [24] in our numerical calculation. In 
this case, the lattice model can preserve the same trans- 
lation symmetry of the discrete version of hik) and turn 
to Eq.(6) in the low energy limit (fc — > 0). We must point 
out that this lattice model to reproduce the low energy 
effective Hamiltonian is not valid for k ^ kp. 

By using the method in Refs. [28, 29], we can obtain 
the surface spectral function. Considering a semi-infinite 
system which has a surface at z = 0, the momentum 
which parallels to the surface fc|| = {kx,ky) is a good 
quantum number, and the partition function of the sys- 
tem with an open surface at z = can be written as: 

Z = f Vi'^V^Pei^p W ^* E £ b^niidt ~ i?fc|, )^„ 

-' [ -^ fc|l n=0 



+ {i^iVk^^^pn+l + h.C.)]} 



where ipn is the wave function for the nth layer, Hk . is 
the intralayer Hamiltonian, V^, , is the interlayer coupling, 
and h.c. means Hermitian conjugate. The recursive inte- 
gration layer by layer gives the following Green's function 
for the surface state: 

G-i(fe||,c.) = Go-i(fc||,c.)-Fj|G-i(fc||,w)\4„, (9) 

where Go(A;|| , uj) = (uj — Hk^.) is the free Green's func- 
tion without interlayer coupling. The Green's function of 
surface state is calculated by the quick iterative scheme 
[30] for T-matrix. Finally, the surface spectral function 
is given in the following form. 



A{ku,uj) = ImTrG(fc||,a;). 

TT 



(10) 



One can also calculate the SDOS by integrating A(fc| | , u) 
over momentum, 



Ps(w) 



d^k 



{2ny 



A{k\\,i 



(11) 



(8) 



In order to agree with the moment independent pairing, 
in the calculation, we make the dimensionless parameter 
A be 0.15, the maximum gap size be about 0.05eV, and 
the truncation of momentum \kx\ = \ky\ be 0.6eV, be- 
cause, as pointed above, the lattice model is not valid for 
k > kp- 

Now we consider the surface spectral function for some 
special pairing symmetries. In the first case, the pairing 
potential is given in Eq.(6), which looks like a direct sum 
of two pairing potentials to describe the ^He-B phase. 
Because, in topological insulator, strong spin-orbit cou- 
pling between different orbits makes the pairing poten- 
tials between different energy bands be complicated, we 
must calculate the topological invariant carefully. We 
can identify the topological invariant by the criterion 
for topological odd-parity superconductor [25-27] which 
shows that the pairing symmetry of Eq.(6) is topological 
non-trivial and we can also calculate the winding number 
directly [15], 

N.^ ^ ^ J d^ke'^''TY[Qid,QkQidjQkQidkQk], 

(12) 
where Qk = 2Pk - 1, Pfe = J^neocc \unik)) {un{k)\ is the 
projector onto the occupied Bloch states. We use the 
later method and find that the winding number is to- 
tally determined by the topology of Fermi surface and 
A'^tu = —1 Sgn(A) if the pairing potential (6) is non- 
vanishing only for a thin spherical shell around the Fermi 
momentum, k ^ kp, which implies that although the pair- 
ing potential is written in four bands (different spins and 
different orbitals) , it can be continuously deformed to the 
weak pairing limit [19] on the Fermi surface. The spec- 
tral function for the Hamiltonian with pairing potential 
(6) (as shown in Fig.l) shows that the bulk state is fully 




-0.8 -0.4 0.0 0.4 



) 


k,/[l/a] 


0.4 


■ 


0.2 


■ 


0.0 


• 


0.2 


■ 


0.4 


■ 



■0.8 -0.4 0.0 0.4 0.8 
k„ [1/a] 




-0.4 -0.2 0.0 0.2 0.4 
k [1/a] 



0.00 



-0.10 -0.05 0.00 0.05 0.10 
CO [eV] 



^^'o.2 M 




-0.8 -0.4 0.0 0.4 
k„[l/a] 



0.8 -0.4 0.0 0.4 O.f 
k„[l/a] 




-0.4 -0.2 0.0 0.2 0.4 
k [1/a] 



0.00 L 



-0.10 -0.05 0.00 0.05 0.10 
(fl[eV] 



FIG. 1: Model calculation of bulk (a) and surface (b) spectral 
function A{k, lj) for BdG Hamiltonian with pairing potential 
Ai,s = Adisig{A,—A,B). (c) Surface spectral function as a 
function of momentum for a; = 0. (d) Bulk (black dash line) 
and surface (red solid line) density of state(DOS). The false 
color mappings of A{k,ijj) in (a), (b) and (c) are in arbitrary 
units. Parameters for model calculation have been given in 
the context. 



FIG. 2: Model calculation of bulk (a) and surface (b) spectral 
function A{k,Lj) for BdG Hamiltonian with pairing potential 
Ai.s = Adia.g{A,0,B). (c) Surface spectral function as a 
function of momentum for a; = 0. (d) Bulk (black dash line) 
and surface (red solid line) DOS. The false color mappings of 
A{k,u:)) in (a), (b) and (c) are in arbitrary units. Parameters 
for model calculation have been given in the context. 



gapped (as shown in Fig. 1(a)) and there is an Andreev 
bound state on the surface (as shown in Fig. 1(b)). Simi- 
lar to the momentum-independent odd-parity pairing po- 
tential A2 in Ref. [25], this pairing potential is inconsis- 
tent with experimental results and there is a minimal of 
the SDOS at w = (Fig. 1(d)), which is not observed 
at the zero-bias conductance peak. One of the explana- 
tion indicates that the bulk band structure is topological 
non-trivial but with some point nodes, which will induce 
a flat dispersion of surface helical Major ana fermions and 
contribute a non- vanishing peak of SDOS at zero energy. 
Among all sorts of momentum-dependent pairing symme- 
tries, there exist some species which possess this property, 
resembling the A4 in Ref. [2^)]. 

In the second case, we consider the pairing potential 
given by Eq.(4) with Ai,s = A diag(A, 0, B) and A"''^^ := 
as other example, and find that the bulk bands have 
two point nodes at fe = (0, i^/^^ — M^/A, 0), and the 
Q-matrix for this pairing potential is well defined in the 
Brillouin zone excluded these two points. In the weak 
pairing limit, the Q-matrix can be written as 



gfc = --Sgn(A)[cos(6')cr^ 



+ sm{e)aa,](E)(Ta-T^), (13) 



near the two singularity points on the Fermi surface, 
where 6 = aTctaTi[Akx/{Bkz)]. Here, we have made a uni- 
tary transformation to express Qk in the eigenvalue rep- 



resentation of h{k) and the a^^z and tq^z in the same form 
as before but with different meanings. Eq.(13) shows that 
there are two ABS on the boundary of xz-plane, which 
is similar to the chiral p-wavc superconductor but time 
reversal symmetry is unbroken here, and these ABS are 
stable if the two point nodes are disconnected. The spec- 
tral function of this pairing potential is given in Fig. 2. 
The SDOS has a non- vanishing value at w = (Fig. 2(d)). 

However, this pairing symmetry seems to be also un- 
likelihood, because there are only two point nodes of the 
bulk bands which breaks the symmetry of D^h group. As 
shown in Fig.2(c), the spectral function of surface state 
for oj = is not invariant under the C3 rotation opera- 
tion in kxky-pleme. As we know, the A4 suggested for the 
pairing symmetry of Cu2;Bi2Se3 in Ref. [24] also exists 
such a deficiency. 

In order to solve this question, in the third case, we 
ask to construct a pairing potential which must satisfy 
the following conditions, (1) it is topological non-trivial, 
(2) its band structure has some point nodes in the k^ky- 
plane and (3) its spectral function preserves the C3 ro- 
tation symmetry, and find that the high order terms of 
momentum arc indispensable. 

Enlightened by the hexagonal warping effects [31] of 
the surface state of topological insulator, we construct 
the following pairing potential, 

Ai(fe) = A[Bfc, + A^3(fc3 _ 2,k.,kl)]ax ® To, (14) 
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FIG. 3: Model calculation of bulk (a) and surface (b) spectral 
function A{k, lu) for BdG Hamiltonian with pairing potential 
given in Eq.(14). (c) Surface spectral function as a function of 
momentum for a; = 0. (d) Bulk (black dash line) and surface 
(red solid line) DOS. The false color mappings of A{k,uj) in 
(a), (b) and (c) are in arbitrary units. Parameters for model 
calculation have been given in the context. 



is bulk gapless when Aj"* ^ and others are zero, and 
the Hamiltonian with only Aj"' ^ can be bulk gapped 
and topological trivial, its winding number N^ = 0. 

In summary, we calculate the spectral function and 
SDOS for three typical momentum-dependent pairing 
potential in a topological insulator, which may act as 
the candidate for the pairing symmetry of supercon- 
ductor Cua-Bia-Scs, we find that similar to momentum- 
independent A2 and A4 proposed in Ref. [25], the pairing 
potentials of the momentum-dependent also permit ABS 
induced by topological non-trivial fully-gapped or node- 
contacted bulk bands, as shown in the first and second 
cases. We point out that the previous topological non- 
trivial node-contacted pairing potentials do not preserve 
the C3 rotation symmetry of lattice structure, and we 
find a solution for this inconsistence, in the third exam- 
ple. 

We must clarify that the node-contacted bulk band 
structure is not the unique explanation for zero-bias con- 
ductance peak, a recent paper [32] shows that a fully 
gapped bulk state with a twisted dispersion of ABS is 
also possible. We hope these different pairing symme- 
tries be judged by future experiment and be helpful for 
the study of the pairing mechanism of Cua;Bi2Sc3. 

Acknowledgement: This work is supported by NSFC 
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where A is a parameter and A = 2eV^^ in calculation, 
for simplicity, we also choose Aj^^^^^ = A2(fe) = 0, as an- 
other example. The term proportional to Bk^ is applied 
to open a gap at the F-point, which can be replaced by 
the A3 in Ref. [25] or some other topological non-trivial 
terms. All of them can give the similar spectral func- 
tions. In addition, we must emphasize that the k^ terms 
are not high order corrections and they are as impor- 
tant as the linear order terms of momentum for the weak 
pairing limit. The spectral function and SDOS are given 
in Fig. (3), now C3 rotation symmetry is preserved (as 
shown in Fig. 3(c)), the net effect of flat helical Majorana 
fermions induced by six point nodes of bulk bands ac- 
cumulates a sharp surface spectral function peak around 
the F-point for w = 0, and this effect is also manifested 
in the SDOS (as shown in Fig.3(d)). 

Finally, we discuss the other choice in Eqs.(3)-(5). In 
general case, the pairing symmetry A^-Jj^, of the anti- 
symmetric orbits is similar to the pairing symmetry A"-'^ 
of the symmetric orbits, we can construct parallel theo- 
ries for pairing potentials which are similar to above ex- 
amples or their combination. For the fully bulk-gapped 
systems, we find that they can deform to each other con- 
tinuously. The difference between A"''^^ and A"-'^ can not 
be distinguished by the shape of spectral function. In ad- 
dition, after calculating the spectral functions for A2(fe) 
at the linear order of momentum, we find that the system 
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